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Abstract 

Geostrophic jets propagating as Alfven waves are shown to arise in a rapidly rotating 
spherical shell permeated by a magnetic field among the transient motions set up by 
an impulsive rotation of the inner core. These axially invariant motions evolve on a 
time-scale which is short compared to the magnetic diffusion time. The numerical 
study is taken as illustrative of a more general point: on such a fast time-scale the 
dimensionless number appropriate to compare the rotation and magnetic forces is 
independent of the magnetic diffusivity in contrast with the often used Elsasser 
number. Extension of the analysis to non-axisymmetrical motions is supported by 
published studies of dynamo models and magnetic instabilities. 
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1 Introduction 



For the last ten years, numerical simulations of the dynamo process in the 
Earth's core have much changed the views on the interplay between magnetic 
and rotation forces. In particular, columnar flows almost invariant in the direc- 
tion parallel to the rotation axis and localized outside the imaginary cylinder 
tangent to the inner core have been found very often even though the Elsasser 
number A, classically used to es timate the ratio of magnetic to rotation force s, 



is of order 1 or larger (see e.g. lOlson et all ( 119991 ): iGrote and Bussd ( I200ll )). 



Alignment parallel to the rotation axis is caused by the predominance of ro- 
tation forces. Accordingly, columnar flows had bee n contemplat ed previously 



in the context of weak- field models (A< 1) alone (IBussd , Il975l ). The "strong 



field regime" (A = 0(1)) was illustrated by mean-field dynamo solutions, in 
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which the azimutal angular velocity showed instead large shears in the di- 



rection parallel to the rotation axis (iBraginskyl . Il978l : iHollerbach and Jones! . 



19931 : iJaultl . Il995h . These early solutions were either steady or slowly varying 
on the magnetic diffusion time in sharp contrast with the current generation 
of dynamo solutions. 

With large magnetic fields, as measured by A, only the geostrophic part of the 
velocity field, symmetric about the rotation axi s, was exp e cted t o be invari- 
ant in the direction parallel to the rotation axis. Braginsky (Il970l ) singled out 



these motions in the context of magnetostrophic equilibrium, characterized by 
the insignificance of inertial and viscous forces compared to magnetic, rota- 
tion and pressure forces. He found that as these azimutal velocities shear the 
magnetic field, they are subject to a restoring force, provided by the magnetic 
field, that ensures wave propagation. This is the mechanism of Alfven waves 
and indeed geostrophic velocities in a rotating spherical shell permeated by a 
mag netic field obey an equation of Alfven wave type save for geometrical fac- 



tors. 



Braginsky! (Il970l ) assigned these torsional Alfven waves to perturbations 
with respect to a slowly evolving basic state characterized by the cancella- 
tion of the total action of magnetic forces on the geostrophic cylinders. That 
description sets the geostrophic velocities apart. My aim, in this paper, is to 
defend another explanation for the emergence of torsional Alfven waves that 
can be generalized to nonaxisymmetric motions, such as the almost axially in- 
variant vortices found in recent geodynamo solutions characterized by strong 
but rapidly fluctuating magnetic fields. Other examples are outlined in the 
discussion part. 

In the next section, I introduce the two dimensionless numbers A and A that 
measure the relative strength of the magnetic and rotation forces, within a 
rapidly rotating body permeated by a magnetic field. I argue that on fast 
diffusionless time-scales, the appropriate number is A. This is illustrated in the 
third section, which constitutes the main body of the article. The competition 
between magnetic and rotation forces is studied in a rapidly rotating spherical 
shell immersed in a magnetic field. Specifically, the axisymmetrical transient 
motions set-up by an impulsive rotation of the inner core are investigated for 
different values of A and the Elsasser number. This is followed by a general 
discussion, where different problems are listed for which A rather than A is 
appropriate to compare magnetic and rotation forces. The paper ends with 
concluding remarks. 



2 Lehnert versus Elsasser numbers 



Elsasser! dl946l ) argued that the magnetic field in the Earth's core saturates 



when the magnetic force becomes comparable to the Coriolis force and sug- 
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where Q is the angular velocity, p is the density and a is the electrical con- 
ductivity. The Elsasser number, 



has subsequently been used to measure the relative strength of Coriolis and 
magnetic forces. In order to derive the relationship (1), the electrical current 
density j is estimated as aUB. This is obviously not valid when magnetic 
diffusion is negligible compared to induction (j <C all B) . 



Conversely, magnetic diffusion does not ent er the p h ysics of plane magneto- 
hydrodynamic waves, of length-scale I that iLehnertl (|l954 ) studied. He used 
another dimensionless number Xo = A -1 , with 



A 



B 
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(3) 



to measure the relative strength of magnetic and rotation forces. The param- 
eter A, hereinafter referred to as the Lehnert number, can be defined as the 
ratio, in a rapidly rotating and electrically conducting fluid permeated by a 
magnetic field, of the period of the inertial waves to the period of the Alfven 
waves. Thus, the typical frequency of diffusionless Alfven waves is Xfl and the 
relationship 

A < 1 (4) 



states that rotation forces dominate over magnetic forces on fast time-scales. 
Cardin et al.l ( 120021 ) argued that both A and A are important to characterize 
geodynamo models. 



In the spherical case, it is convenient to specify A using the outer radius a as 
the length-scale I in the definition (3). Denoting by E M the magnetic Ekman 
number r]/Qa 2 and by E the ordinary Ekman number z//f2a 2 ,the relationship, 



E M + E < A, 



(5) 



ensures that Alfven waves are not rapidly damped by either magnetic or vis- 
cous diffusion [y and rj are respectively the viscous and magnetic diffusivities). 



3 



Figure 1. Field lines of the imposed poloidal axisymmetric field. 



Indeed, the ratio of A to E + E M is the Lundquist number S which indi- 
cates how far A lfven waves propagate before they are quenched by diffusion 
( IRobertd . 119671 ) . The two numbers A and A are related through the magnetic 
Ekman number: 



A = AE 



M- 



(6) 



The value of A appropriate to the Earth's core is of the order of 3. x 1CT 5 — 
2. x 10~ 4 . Indeed, A = 1CT 4 corresponds to a magnetic field strength in the 
core interior of the order of 3. mT. Usua lly quoted values range from 1. mT 
( IChristensen and Aubertl . l2006l ) to 4. mT f jStarchenko and Jonesl . l2002l ). Using 
E M = 4. x 10" 9 yields S of the order of 10 4 - 5. x 10 4 . 



3 Axisymmetric motions spawned in a spherical cavity by a sudden 
impulse of the spin of the inner core 



The change of the relative strengths of the magnetic and rotation forces ac- 
cording to frequency is well illustrated by the contrast between transient and 
steady flows in a differentially rotating spherical shell in the presence of a mag- 
netic field. Static solutions have been published in the case of dipolar magnetic 
field and small differential rotation, for which the structure of the flow has been 



descr i bed according t o the E lsasser number f jHollerbachl . Il994l : iDormy et al. 



1998 ). Kleeorin et al. (Il997l ) have theoretically investigated steady linear so- 



lutions when the imposed magnetic field is potential and has dipole parity. 
They have identified several asymptotic regimes according to values of the El- 
sasser number in the small Ekman number limit. Let us now study transient 
structures. 
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3. 1 Model and governing equations 



Consider an electrically conducting homogeneous fluid occupying a spherical 
shell that is immersed in an imposed steady magnetic field. The ratio of the 
inner shell radius b to the outer radius a is b/a = 0.35 as in the Earth's core. 
The solid inner core has the same electrical conductivity as the fluid and the 
outer boundary is insulating. The fluid is rotating with the constant angular 
velocity fl. The imposed magnetic field is chosen as: 



B = B Vx (Aet) (7) 
A=(j 1 (Jhir) - 0.3j 1 (A 2 r))P 1 1 (cos^) - 0.2 j3 (f3 31 r)P^cos6) (8) 



where (r, 9, 0) are spherical coordinates, are Legendre functions, is 
the set of spherical Bessel functions of the first kind and (3i n is the nth root 
°f ji-i(P) — 0. The basic field, which is shown in figure 1, has been chosen 
with the aim of modelling torsional oscillations in the Earth's fluid core. In 



this context, the important quantity (|Braginskyi . Il970l ) is 




MK«) = -T f Bidzd<P , (9) 



evaluated on geostrophic cylinders of radius s and of top and bottom z 
coordinates respectively zt and zb- Obviously, the choice (8) is arbitrary. 



In con trast with the Earth's case and with the basic state used by iBraginsky 



(119801 ) to model torsional oscillations, {B^ } vanishes at s = a because of the 
imposed dipole symmetry with respect to the equatorial plane. Axisymmetry 
makes {B^} = at s = in contrast with the geophysical case again. In 
view of the present study, the main characteristics of the field B defined by 
(8) are that it is neither parallel to the rotation axis nor rapidly decreasing 
with radius as are current-free dipole fields. The results presented below do 
not depend on the details of the geometry of B. 

Study the evolution of the velocity field u and of the magnetic field deviation 
b after an impulsive increase of the angular rotation ojb of the inner core, 
postulating symmetry about the axis of rotation and dipole symmetry about 
the equatorial plane: 



u r (r, 7r - 9, </>) = u r (r, 9, 0), u e (r, ir - 9, 4>) = -u e (r, 9, 0), 

u^r, 7i - 9,<fi) =u (f) (r, 9,<j)), b r (r,Tc - 9,<fi) = -b r (r,9,<p), (10) 

b e (r, 7i - 9, <p) = b e (r, 9, (ft), b^r, tx - 9, 0) = -b^r, 9, 0). 
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The amplitude of the initial impulse Q& is assumed to be small enough so that 
the subsequent evolution of the dynamics is independent of within a scaling 
factor. Using B as unit of magnetic field, a as length-scale and a(/xp) 1/ ' 2 / B 
as unit of time, the fields u and b are governed by the following linearised 
equations in the fluid region: 



— + 2\- 1 e z x u = - Vp + (V x B) x b + (V x b) x B 

at 

+P m AA- 1 V 2 u, (11) 
^ = Vx(uxB) + AA-Vb, (12) 

where P m = v/rj is the magnetic Prandtl number. The field b is defined also 
in the inner solid region where: 

f = AA^V 2 b. (13) 



Note that the steady-state solutions depend only on the two parameters A 
and E. The velocity boundary conditions 

u = 0, r = a (14) 
u = s u) b {t) = s Qb S(t — t )> r = b (15) 

are written using the Dirac 5 function and are appropriate to rigid boundaries. 

The set of equations (11) and (12) is discretized and time-stepped from an 
initial state of rest: 

u = b = (16) 



The Dirac S function is approached as: 



1 



e -(t-to) 2 A ( 17 ) 



It is a result of the simulations below that the solutions are independent of 
the parameter e provided that its value is set small enough. 

A poloidal/toroidal decomposition 



u = u (j) e ( f > + V x (w p e ) (18) 
b = 6^ + V x (bpe^) (19) 
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(a) 



(b) 




(c) 



(d) 



Figure 2. Contours of the induced azimutal magnetic field for A = 1.72 x 10 , 
A = 0.52, E = E M = 5.7 x 10~ 8 drawn at t = 8.6 x 10~ 3 (a), t = 1.72 x 10~ 2 (b), 
t = 3.44 x 10~ 2 (c) and t = 6.88 x 10 _2 (d). The contour intervals and the frame size 
are identical in all frames. The inner sphere boundary is indicated with a thick line. 

is employed. The variables are expanded in associated Legendre functions, i.e. 



and then discretized in radius. The minimum truncation level Imax is 120 
whereas at least 450 unevenly spaced points are used in the radial direction. 

3.2 Formation and propagation of geostrophic jets 

Let us examine a typical sequence of solutions for a small value of A. Following 
the initial impulse, an almost geostrophic shear sets up, after a few revolutions, 
at the cylindrical surface tangent to the inner core, hereafter referred to as tan- 
gent cylinder. Induction of an azimutal magnetic field localized at the tangent 
cylinder (compare the snapshot (a) to the snapshot (b) in figure 2 occurs about 
the end of this period during which the velocity field becomes axially invariant. 
It starts from a source at the equator of the inner core. The last two panels 



Imax 



= E ul 4>( s ) P 2i+i(cos9) 



(20) 
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(a) 



(b) 




(c) 



(d) 



Figure 3. Contours of constant angular velocity for A = 1.72 x 10 4 , A = 0.52, 
E = E M = 5.7 x 10~ 8 drawn at t = 8.6 x 10~ 2 (a), t = 0.26(b), t = 0.52(c) and 
t = 1.03(d). The contour intervals are respectively Aw, Aoj/2, Au/5 and Aw/10 in 
the frames a, b, c and d in order to compensate for the attenuation of the velocities. 

of figure 2 illustrates the following period during which meridional electrical 
currents parallel to the tangent cylinder intensify and loop further and further 
away from the tangent cylinder. This is the most noticeable feature before 
the geostrophic shear splits up. The outer shear readily transforms into a jet 
propagating away from the tangent cylinder towards larger cylindrical radii. 
Thereafter, the velocity within the jet becomes more and more invariant along 
z as time elapses (figure 3). On the other side of the tangent cylinder, a second 
shear propagates towards the axis of rotation. Its propagation velocity slows 
down as B s decreases to on the axis. In the event, the inner shear transforms 
also into a jet. The comparison (figure 4) with a second sequence of solutions 
for the same value of A but for A multiplied by a factor of 12.5 indicates that 
the dynamics outside the tangent cylinder is almost independent of A. The 
flow remains geostrophic even though A is 0(1). Note that steady flows do 
not reproduce this feature. Figure 5 shows zonal flows driven by rotating the 
inner sphere at a constant rate for the two values of A used to calculate the 
transient solutions. For the largest value of A, the angular velocity contours 
are not parallel to the rotation axis and tend instead to follow the magnetic 
field lines, as prescribed by Ferraro's law of isorotation. Steady solutions are 
established after a period lasting a few time units A -1 !] -1 during which the 
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(a) (b) 

Figure 5. Steady azimuthal flow induced by rotating the inner and outer boundaries 
at slightly different rates. Contours of constant angular velocity for E = 4.5 x 10~ 7 
and A = 0.52 (a) A = 6.5 (b). 

flow is geostrophic. 

The outer geostrophic jet has finite width 8 in the limit e — > (see expression 
(17) for the definition of e). Investigating the variation of S with A, E, and 
Em gives an useful insight into the mechanism of generation of the geostrophic 
motions. Here S is arbitrarily defined as the distance along s between the two 
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Figure 6. Scaling of the width S of the outer geostrophic jet with respect to the 
number A. The magnetic Prandtl number is P m = 1 and the Ekman number is 
E = 5.7 x 10~ 8 (crosses) or E = 2.85 x 10 -8 (+ signs). A line of slope —1/4 is shown 
for comparison. 
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Figure 7. Scaling of the width 5 of the outer geostrophic jet with respect to the 
Lundquist number = X/Em- The magnetic Prandtl number is P m = 1 and A is 
4.3 x 1(T 5 (O), 8.6 x 10" 5 (x), 1.72 x 1(T 4 (o), 3.44 x 10" 4 (+), 6.88 x 1(T 4 (□). 
A line of slope —1/4 is shown for comparison. 

cylinders, on both sides of the geostrophic jet, where the angular velocity has 
half its maximum value. For the range of parameters that has been extensively 
explored, the outer layer is always well characterized from t = 0.17 onwards. 
Results are reported for this time. Keeping E constant and P m = 1, it is found 
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Figure 8. Thickness S of the outer geostrophic jet with respect to P m l for P m < 1. 
and various A: A = 4.3 x 1(T 5 (circles), A = 8.6 x 1(T 5 (crosses), A = 1.72 x 1(T 4 
(+ signs), A = 3.44 x 1CP 4 (squares). 



that S varies as A~ 1//4 (see figure 6). For fixed values of A and P m 
as E^ 4 . This is illustrated by the figure 7. Indeed, assuming 5 
power law 5 ~ E 1 ^ can simply be written: 



(*) 



-1/4 



1, 5 varies 
the 



(21) 



using the Lundquist number = X(Em)' 1 - Thus, the width 5 is independent 
of the angular velocity Q. We are interested by results for P m < 1 since v <C rj 
in the geophysical case and in laboratory experiments as well. Decreasing P m 
from P m = 1, a slight dependence of S on P m is found (figure 8). Extension of 
the relationship (21) to solutions for P m < 1 is supported by these results. 

The outer geostrophic jet is radiated from the tangent cylinder, which touches 
the inner core on its equatorial circle. There, both the rotation vector and the 
magnetic field are parallel to the inner core surface and the Ekman-Hartmann 
viscous boundary layer adjacent to the inner core is singular. Thus, it is in- 
structive to investigate the influence of the strength of the magnetic field at 
the singularity. For this purpose, an axial uniform field can be added to the 
magnetic field defined by (8). It is found that the jet is much thickened if the 
two fields cancel out at the singularity. Then, the relationship (21) does not 
hold. On the other hand, it is also found - in the narrow parameter range 
which has been investigated - that the layer shrinks as [B z (b, O)]" 1 ^ 4 as the ax- 
ial magnetic field adjacent to the equatorial ring of the inner core is increased. 
Taking this result at its face value, the magnetic field strength entering the 
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relationship (21) would be B z (b,Q). Putting these results together, it appears 
that the magnetic structure adjacent to the equator of the inner core plays 
a significant role in the emergence of the two propagating shear layers. The 
transformation of the outer shear layer into an independent jet detached from 
the inner core is promoted by the axial magnetic field. 

Once the outer j et is formed, its time evolution is given by the equations of 



Braginskyl (119701 ). which satisfy angular momentum conservation. It is possible 



to estimate the angular momentum A(t) = s 2 Vl — s 2 5u ( f,(s) carried by the 
outer jet from the width S. For the solutions that have been investigated, A 
does not change throughout the propagation of the jet. The strength B s of the 
magnetic field sheared by the jet decreases with s to at s — 1. Thus, the jet 
slows down as it approaches the outer sphere equator, which it never reachs. 

The scaling (21) has implications for the coupling betwe en the axial rotation 



of th e Earth's solid inner core and torsional oscillations (jBuffett and Moundl . 



20051 ). The time unit A 1 Q 1 corresponds to 1-10 years for geophysical appli- 
cations. A characteristic time (S 1 ^ (Afi)" 1 can be derived from the length 
S and the Alfven wave velocity Xfla. It corresponds to the inner core rotation 
period below which dissipative processes at the tangent cylinder are impor- 
tant. Using a geophysical estimate for , the coupling mechanism presented 
above between the rotation of the inner core and torsional Alfven waves is 
found to be efficient on periods longer than a few months. Investigation of 
magnetic fields with non-dipole symmetry will be a natural follow-up of this 
study. 

Finally, keeping constant and increasing A, it is found that the structures 
radiated from the tangent cylinder lose their geostrophic character for A ~ 

io- 2 . 



4 Discussion 



Let us now examine to what extent the parameter A is also appropriate to 
comment on the geometry of fast motions in other problems characterized by 
rapid rotation and magnetic field. 



4-1 Axially invariant hydromagnetic instabilities occurring at small Lehnert 
number 



We can discriminate between two approaches that have been followed to study 
the stability of a magnetic field in a rotating sphere according to the parameter, 
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either A or A, used to compare magne tic and rotation forces. iMalkud fl 19671 ) 
recently followed by lZhang et al.1 (120031 ) studied hydromagnetic waves in a non- 
dissipative fluid (A — > oo). As a result, they wrote the condition for stability 
as a relationship i nvolving the Lehnert nu mber A, which has to exceed values 
of the order unity. IZhang and Fearnl (119941 ) focused on the rapid rotation limit 
(A — ► 0) instead. Then, the onset of instability occurs for a critical value of the 
Elsasser number A c . In accordance with the above discussion on the geometry 
of the motions in the limit (A — > 0), they found that the (non-axisymmetric) 
instability is characterized by nearly two-dimensional columnar fluid motions 
despite A c being O(10). 



This result stands when the instabilities are thermally driven. IZhangi (119951 ) 
focused his study of rotating convection in the presence of an axisymmetrical 
toroidal magnetic field on the limit (E — > 0), which amounts to (A — > 0) for 
finite values of the Elsasser number A. The magnetic Prandtl number is set to 
1 and the control parameters are thus A together with a Rayleigh number. The 
fluid motions, in the limit (E — > 0), are almost two-dimensional showing only 
slight variations along the direction of the rot ation axis whilst A is 0(1 — 10). 
The solutions of lOlson and Glatzmaien (I1995h (ou t side t he cylindrical surface 



tangent to the inner core), IWalker and Barenghil (|1997) (th eir figures 6f, 7f, 
8f, 9f) for different basic states and Zhang and Jones! T 19961 ) (their figure 4) 
for n/r] <C 1, where k is thermal diffusivity, all present similar features. 



4-2 Columnar flow structure in geodynamo models 



Obviously, the numbers A and A are less directly relevant to studies of dynamo 
simulations than to investigations of models with imposed large-scale magnetic 
field. These two estimates of the magnetic field strength come out as output 
of the numerical runs instead of being among the initial parameters. It is nev- 
ertheless true that realistic values of A are reached in numerical models of the 
geody namo as A does not involve diffusivities. Thus, IChristensen and Aubert 
( 120061 ) conducted a statistical analysis of a set of geodynamo models and esti- 
mated a parameter defined as A, using the shell depth as the length-scale I in 
(3). Their results correspond t o A varying from 7 x 10~ 3 to 3 x 10~ 2 , keeping 
the core radius as length-scale. IChristensen and Aubertl (j2006l ) found that the 
narrow range of A values contrasts with the wide variations of A. They also 
remarked that measuring the relative strength of magnetic and rotation forces 
acting in geodynamo models with the Elsasser number A does not reflect the 
fact that the Lorentz force depends on the length scale of the magnetic field 
whereas the Coriolis force is independent of the length scale of the velocity 
field. Conversely, comparing typical periods of the Alfven waves to typical 
periods of the inertial waves shows that the relative importance of the mag- 
netic force augments with decreasing length scales (see the expression (3) of 
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A ) . In their previous systematical parameter study, IChristensen et al 



(1393 ) 

had found only one case with strong deviations from columnarity (A = 14, 
E M = 6 x 10~ 5 , A = 3 x 10~ 2 ). Together, these results are consistent with the 
statement that the extent to which rotation affects the structure of the flow 
depends on A, the motions remaining columnar up to A = 0(3 x 10 -2 ). 



4-3 Torsional oscillations in convective dynamo models 



Convection columns can excite time-dependent geostrophic motions in dy- 
namo mo dels through magnetic and Rey nolds stresses. This has been illus- 
trated by iDumberry and Bloxhaml (120031 ) . They separated the axisymmetric 
azimu tal velocity field obtained from the geodynamo model of lKuang and Bloxham 
( 119991 ) - E = E M = 4 x 10~ 5 , P r = 1 and stress-free boundary conditions - into 
a mean flow plus a fluctuating component. They showed that the quasi-static 
azimutal winds have lar g e gra dients in the z-direction. On the other hand, 
Dumberry and Bloxham (120031 ) emphasized the axial invariance of the time- 



varying zonal flows. Their finding that, outside the tangent cylinder, the whole 
length of the geostrophic cylinders accelerates azimutally as if they were rigid 
on time-scales r ~ 0.1 r D is in line with the small value of A in this numerical 
experim ent (rn magnetic diffusion t ime). Indeed, using B ~ 2(2fi / up?7) 1 / 2 (see 
fig. 10 of lKuang and Bloxhaml (Il999h ). we infer A ~ 10 and A ~ 2xlO~ 2 .These 
geostrophic motions are not Alfven waves as Reynolds stresses and viscous 
forces are as important as the magnetic forces in the balance of forces acting 
on the geostro phic cylinders. More rec ently and with Earth-like no-slip bound- 
ary conditions. iTakahashi et al.l (120051 ) argued indeed that the Ekman number 
has to be decreased down to E = 8 x 10" 6 to make the viscous torque acting 
on the geostrophic cylinders negli gible and the magnetic torque predominant. 
Finally, for the same value of E as ITakahashi et al.l (120051 ). but with stress-free 



bound ary conditions and small Prandtl number P r = 0.1, iBusse and Simitev 
(120051 ) found a dynamical state where the magnetic torques on geostrophic 
cylinders account for most of the geostrophic accel eration. Extracting t he av - 
erage magnetic field strength from the figure 18 of IBusse and Simitevl (120051 ) 
gives A ~ 5 x 10~ 3 - well in the domain A < 1 - and A ~ 3. The result 
that magnetic forces dominate over Reynolds stresses in the balance of force 
acting on the geostrophic cylinders can be related to the observation that the 
magnetic energy is much stronger than the kinetic energy in this solution. 
Thus, sequences where geostrophic motions behave as torsional oscillations 
begin to be detected in convective dynamo models. That requires A < 1 - 
to obtain time-dependen t geost rophic motions, observed for A ~ 2 x 10~ 2 by 



Dumberry and Bloxhaml (120031 ) -, small E and, presumably, kinetic energy 



weaker than magnetic energy. In these studies, there is no evidence that the 
quantity {-B 2 } measuring the intensity of the magnetic field sheared by the 
geostrophic motions and the geostrophic velocities evolve on separate time- 
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scales. Further work is needed to decide what kind of models (fully consistent 
dynamo models with poor separation of scales versus models incorporating a 
steady field) better describes the Earth's core dynamics. 



4-4 Torsional oscillations and Taylor states 



In this article, torsional oscillations are considered as part of the rapid mo- 
tions that a re do minated by rotation because magnetic diffusion is negligible. 



Braginskyl (119701 ) had a different line of arguments. He attributed torsional 



oscillations to departures from a dynamic equilibrium w here the net torque 
of the Lorentz force on any geostrophic cylinder is zero (iTaylorl . 1 19631 ) . This 
condition has to be met, in spherical geometry, when only the Coriolis, pres- 
sure, buoyancy and magnetic forces are taken into account (MAC balance). 
It is frequently referred to as a "Taylor state". It describes a dynamo regime 
on the long time-scale for which magnetic diffusion is important. Reinstating 
the acceleration of geostrophic motions du$(s)/dt in the equation for azimutal 
velocities, i t has been possible to exhibit invis cid solutions of the model-Z of 



Braginskyl ( 119781 ) that are in a Taylor state f jJaultl . Il995l ). However, this is 



almost the unique instance where a connection between torsional oscillations 
and idealized Taylor states has been vindicated. Geodynamo numerical mod- 
els showing torsional oscillations that keep bringing back the magnetic field 
towards a Taylor state have not yet been found. The two viewpoints differ in- 
sofar nonzonal rapid motions are considered. I envision here that they are also 
constrained by rotation being almost ^-invariant whereas iBraginskyl (119701 ) 
made no predictions on the geometry of these motions. 



5 Concluding remarks 



Focusing a numerical study on the transient motions spawned by an impulse 
in the rotation of the inner boundary of a rapidly rotating spherical shell 
immersed in a magnetic field with dipole symmetry, I have documented the 
emergence of geostrophic jets from the cylindrical surface that touches the 
inner core at its equator, irrespectively of the value of the Elsasser number. 
Both the poloidal motions, of which the vorticity is aligned along e^, and the 
toroidal motions with shear in the z direction are rapidly eliminated. The 
geostrophic layers travel with the velocity Xa fl and are govern ed, outside the 
tangent cylinder, by the equations written by IBraginskyl ( 11970l ). The jet width 

scales as (S^j . This estimate yields the frequency below which oscillations 
of the solid core are communicated to torsional Alfven waves in the fluid shell. 
Using Earth-like parameters, it corresponds to a period of a few months. This 
study gives an illustration of the key role played - for fast flows - by the 
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param eter A - independent of magnetic diffusivity - put forward by iLehnert 



(119541 ) . Conversely, A is not appropriate to study steady solutions as it cannot 
be derived from the two parameters A and E that characterize the static 
problem. 

In the same spirit, I have been able to base a discussion of earlier numerical 
studies of magnetic instabilities and dynamos in rotating shells on that pa- 
rameter A. I suggest that the smallness of A in some of these studies is the 
reason for the occurrence of columnar motions aligned parallel to the axis of 
rotation and also of geostrophic flows evolving as Alfven waves. From these 
earlier studies, I anticipate that the results presented here can be extended to 
the nonaxisymmetric case. 

Thus, the parameter A, instead of the Elsasser number A , is the appropri- 
ate parameter to compare magnetic and rotation forces when flows evolving 
on time-scales much shorter than the magnetic diffusion time are considered. 
As the value of A appropriate to the Earth's fluid core is O(10~ 4 ), I advo- 
cate that motions in the core interior with fast diffusionless time-scales are 
approximately z-independent and columnar with vorticity aligned parallel to 
the rotation axis. That paves the way for dynamical studies of the flows re- 
sponsible for the secular variation of the Earth's magnetic field, generalizing 
to all fast motio ns what has already been achi eved for the geostrophic, axially 



symmetric ones ( IZatman and Bloxhaml . 119971 ) . 
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